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Abstract. 

The advent of quantum computing has heralded a renewed interest in physical 
memories - physically realizable structures that offer reliable data storage with error 
correction only at the point of access. Here, we examine a model of a classical physical 
memory, capable of storing a classical bit, based on the ferromagnetic Ising model on 

1- and 2-dimensional square lattices. We make use of Monte-Carlo simulations as well 
as analytic solutions in the high temperature limit. At high temperatures, both 1- and 

2- D Ising models behave like a non-interacting model, albeit with a reduced effective 
number of spins and a reduced effective spin-flip rate. In agreement with general 
arguments, we confirm numerically that (i) the information storage time is independent 
of system size in ID, (ii) in two dimensions, at temperatures below the Onsager phase 
transition temperature Tc, the information storage time increases exponentially with 
system size. Interestingly, some benefit in increasing system size is found even above 
Tc, although the increase in lifetime is not exponential. 
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1. Introduction 

The advent of quantum computing [Ij has heralded a renewed interest in self-correcting 
physical memories - physically realizable structures that offer reliable data storage with 
error correction only at the point of access. Physical storage of quantum bits presents 
many challenges, and model memory structures such as the surface code Hamiltonian [2] 
have been devised to meet these. In classical computing, by contrast. Nature has pro- 
vided us with a physical memory: the ferromagnet, which is the basis of magnetic data 
storage [3]. This article investigates quantitatively the use of the ferromagnet as a physi- 
cal memory. We hope that the understanding gained of this simple case will also provide 
a starting point for a deeper understanding of quantum memories. 

The Ising model, proposed by Lenz in 1920 [1], considers an anisotropic ferromagnet 
as a D-dimensional lattice of spins, each of which may take values Sj = {1,-1}. The 
energy of this Ising lattice is given by: 

if = - J ^ SiSj - ^ hiSi (1) 

<i,j> i 

where the first sum is over nearest neighbour pairs. J represents the strength of the 
magnetic interaction between neighbouring spins, while h is the external magnetic field. 
For a ferromagnet, J > 0. In this work, we take h = 0, J = 1. For h = 0, H is 
minimised when all spins are equal, so the Ising lattice has two energy minima. 

Initial investigations asked whether or not an infinite Ising lattice displays a phase 
transition. This would correspond to the transition from para- to ferro- magnetism as 
the temperature of a magnetic material is decreased. Ising [5] showed that there is no 
phase transition for D = 1. Kramers and Wannier [6l[7j later showed that a phase 
transition does occur for = 2. Onsager [8] identified the critical temperature at which 
this transition occurs as: 

9 T 1 kT 

tanh = ^ ^ ^ ^ 2.269185 (2) 

The critical temperature is equivalent to the Curie temperature of a magnetic material: 
for T < Tc spontaneous magnetization occurs. The magnitude of the spontaneous mag- 
netization was calculated exactly by Yang [S] in 1952. 

The energy landscape of an Ising lattice is the basis of classical magnetic data 
storage. We model the storage process by attempting to store a single bit, x, using an 
Ising lattice of spins with side length L. Each energy minimum corresponds to a 
different value of the bit. For example, we might choose: 

Si = — IVz <^=^ X = 

(3) 

Si = l\/i <^==^ X = 1 

Thermal spin fluctuations may corrupt our lattice. Fortunately, these will carry an 
energy penalty, so the system will tend to correct such fluctuations. In this way, the 
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Ising lattice acts as a self-correcting physical memory. After some time, we attempt to 
retrieve our stored bit by majority voting. This corresponds to a final error correction 
step at the point of readout. Such a step is always necessary when evaluating the per- 
formance of a self-correcting memory based on an error correcting code, as has been 
noted by Pastawski et al. ^U\. For example, if more than hali[§| of the spins have the 
value -1 we infer that x = 0. From an information-theoretic viewpoint, the Ising model 
is a physical realisation of a repetition code. Likewise, the surface code Hamiltonian 
represents a physical realisation of a quantum error correcting code. 

The energy penalty during the evolution of the system occurs on the perimeter of 
regions of different spin. For a ID Ising lattice, a region of incorrect spin may therefore 
grow without further energy penalty. However, growth of incorrect spin regions in 
lattices with D > 1 requires further energy input from the environment. This argument, 
proposed by Peierls [11], leads one to expect that the ID Ising lattice is a poor physical 
memory, whereas lattices with D > 1 are good physical memories at sufficiently low 
temperature. One aim of this paper is to explore quantitatively these ideas via Monte- 
Carlo (MC) simulations. 

An analogous effect is seen in quantum physical memories, in particular in the toric 
codes described in reference [2]. In a 2D toric code, the logical operators that take the 
system from one ground state to another are string-like. This means that, once created, 
quasi-particles representing the propagation of these errors can travel without energy 
penalty. Consequently, the 2D toric code is not a good physical memory. In the 4D 
toric code Hamiltonian, by contrast, logical operators are membrane-like, and therefore 
the propagation of errors is accompanied by an additional energy penalty, and a good 
physical memory is expected to result |2l|12]. Recently, a 3D stabilizer code Hamilto- 
nian with no string-like logical operators has been proposed, which does energetically 
penalise additional errors and therefore may serve as a useful quantum information stor- 
age device p^HTS] . 

We define the fidelity F{t) of a memory as the probability that the correct bit will 
be retrieved at time t, given that the bit was encoded perfectly at t = 0. For an Ising 
lattice, F{t) is expected to depend on dimensionality, geometry and temperature. In 
this article, we will develop a quantitative model for F{t) for Ising lattices with free 
boundaries and no external field. We will use MC simulations to verify the model and 
investigate its implications for ID and 2D lattices at various temperatures. 

We observe two distinct mechanisms of fidelity loss in an Ising lattice, depending 
on its temperature. For very high T, the lattice will decay into a disordered state, 
with no large regions of equal spin. For very low T, the lattice will suffer only small 

§ In the unlikely event that exactly half of the spins have the value -1 we have two options: (1) Choose 
the value of x randomly or (2) Declare that the data has been corrupted. The latter is mathematically 
equivalently to inferring the wrong value of the bit: in either case, the memory has failed. We will use 
whichever of these is mathematically convenient. 



The Ising ferromagnet as a self-correcting physical memory: a Monte-Carlo study 4 

fluctuations from its initial state for a long period of time, until the unlikely event that 
a fluctuation grows to cover the whole of the lattice. The lattice will then remain in 
this new state, with the vast majority of spins incorrect, until another such fluctuation 
reverts the lattice to its original state. For intermediate T, at sufficiently large times 
the lattice is composed of a few regions of equal spin. These ideas are seen also in work 
on the kinetics of an Ising lattice in an opposing magnetic field (i.e. with spins initially 
aligned against the field). This area has benefited from extensive research for a number 
of years, see for example references [T6ti30] . In this environment, droplets of spin aligned 
with the field will form, but due to the energy penalty at the perimeter small droplets 
are likely to shrink. However, if a droplet reaches a critical size it will grow to fill the 
entire lattice. When h 0, the lattice will likely remain in this configuration (the 
minimum energy state) for an extremely long time. In our case, the lattice has a free 
choice of energy minimum, and so will continue to behave stochastically. 

Ising lattice fidelitjjJJ] was studied using MC simulations by Kolesik et al [23] for a 
lattice with D = 2 at temperatures T < Tc in an opposing field with periodic boundary 
conditions. The Ising lattice with free boundaries in an opposing field has also been 
investigated, for example by Cirillo and Lebowitz [25], but the fidelity is not consid- 
ered in this work. Cirillo and Lebowitz demonstrate that the use of free rather than 
periodic boundaries is highly significant to the dynamics of the Ising lattice. With free 
boundary conditions, the critical droplet size, and hence the lifetime of the metastable 
state, is lower. It was also noted that the critical droplet develops in a corner of the 
system, where the energy of the perimeter is lower. Real ferromagnets, although they 
may display an exceedingly large value of A^, have free boundaries. 

In the next section we consider the dynamics in the limiting cases of high and low 
T. We then describe our procedure for simulating the time evolution of the fidelity in 
the case of arbitrary temperatures, and propose a model function for this time evolution, 
inspired by the high-temperature limiting behaviour. In section 3 we present and discuss 
the results of our Monte-Carlo simulations. In section 4 we conclude and point out some 
questions for future study. 

2. Dynamics 

2.1. Limiting Cases 

In general, the evolution of the lattice under the influence of thermal noise will depend 
on the details of the system-bath interaction. In order to make progress, we do not 
assume a specific form of this interaction, but rather model the noise as a stochastic 
process that acts locally to flip spins and yields a long-time solution that is consistent 
with a Boltzmann distribution at temperature T. We first consider the limiting cases 

II This work studies the probability that the sign of the magnetization remains unchanged, as this is 
measurable with magnetic force microscopy experiments. This is equivalent to fidelity as defined here. 



The Ising ferromagnet as a self-correcting physical memory: a Monte-Carlo study 5 

of very high and very low temperature. We begin our investigation by considering 
non-interacting spins. This is the hmit as T — oo of any Ising lattice with spins. 
At t = 0, we encode a; = 1 i.e. Si = IVi. At some later time, the spin Sj has changed 
value Ui times, where rij is modelled by a Poisson distribution with some characteristic 
flipping rate A per unit time. If n is even, Sj = 1; if n is odd, Sj = —1. Therefore: 



The fidelity is the probability that fewer than half of the spins have flipped an odd 
number of time^ which can be obtained using the binomial distribution: 



In the very low temperature case, the lattice is usually found in a state in which 
almost all spins are equal because the energy barrier between these minima is very 
high. As explained above, loss of fidelity is caused by the non-zero probability of the 
lattice transitioning between the minima. Therefore an approximate model for the low 
temperature case considers the whole lattice as a single spin with a very low A. This 
phenomenon is known as superparamagnetism [31] and presents a challenge for mag- 
netic data storage [321431] . This simple interpretation gives an exponential form of the 
fidelity (equivalent to setting A^ = 1 in (5) ), as noted by Kolesik et al [23] . 

2.2. Monte-Carlo Simulations at Arbitrary Temperatures 

We used an MC process to simulate the evolution of the Ising lattice. At each MC 
time-step, we updated a single randomly chosen spin. We flipped the sign of the spin 
with a probability given by Glauber dynamics: 

P{fl^p) = (6) 

1 + e fcT 

where AE = energy of lattice with spin flipped — current energy of lattice. In order 
to make meaningful comparisons between systems of different sizes, physical time is 
modelled as t = _ We approximated the fidelity by repeating each simulation 

M = 10^ times and finding the fraction of these that encoded the correct value of x. In 
the non-interacting case, we have always P{flip) = 0.5, so each spin flips on average 0.5 
times per physical time-step, and we expect A = 0.5. At any given time, each of our M 
runs has F = 1 or F = 0. Therefore F{t) has a binomial distribution with M trials and 
probability of success F{t), normalised by dividing by M. This allows us to find the 
statistical error in F{t). Because M is finite, we have also a discretization error. Errors 

% Here we have chosen option (2) for the case of a 50/50 spht. 




(4) 
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in the fidelity calculated from the MC data were approximated by taking the largest of 
the discretization error and the statistical error: 

1/2^ 



MC simulations of the non-interacting case are an excellent fit to the analytical 
result, as shown in figure 1. As we discuss further in later sections of this paper, the basic 
form of the fidelity curve for an Ising lattice is a deformation of that for non-interacting 
spins. In particular, we note that, even at low temperature and in both dimensionalities 
simulated, the initial region of the fidelity curve with negative curvature (seen clearly in 
the non-interacting case) is still present. Therefore a simple exponential model cannot 
capture all the features of fidelity decay. 

Fidelity of 100 non-interacting spins 
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Figure 1. MC data (blue with vertical error bars) and model (red) for the fidelity of 
100 non-interacting spins, — 1554; Degrees of Freedom = 1550 
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2.3. A Mathematical Model for F{t) 

Observations of the evolving Ising lattice verify that, in general, the lattice exhibits 
several spin regions (see figure 2) , which may move around the lattice and change sign. 
We therefore postulate that the Ising lattice can be approximately modelled using (5) 
by allowing A and the effective number of spins iVeg to vary. iVeg is interpreted as the 
average number of spin regions in the lattice. 
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Figure 2. A square Ising lattice with N = 100 and kT = 2.5 after 177000 MC steps. 
Initially all spins were +1 (red). Our model predicts N^s — 5.37 ± 0.09 for this lattice. 
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To allow continuous variation of N^g, we used a Gaussian approximatioEEl for the 
binomial distribution. Our proposed model, then, is: 



/i 



a 



2^efT(l - 



-2Xt\ 



-2\t 



(8) 



)(l + e--*)] 



-2At\]l/2 



Figure 3 shows how the shape of F{t) varies with N^s and compares the Gaussian 
approximation to the exact binomial form. We can see that decreasing Nes decreases 
the size of the flat region at low t with respect to the length of the decay, and overall 
quickens the decay of the fidelity. For a given lattice, however, a decrease in N^s is 
accompanied by a decrease in A, as the larger spins regions flip less frequently. A 
decrease in A preserves the shape of the curve, but stretches it along the time axis. 



The model with X = 0.1 
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Figure 3. Binomial (lines) and Gaussian (points) form of the model with A = 0.1; 
Nf.ff = 399 (solid line, crosses) and A^e// = H (dashed line, circles) 



3. Results and Discussion 

We simulated ID and square 2D lattices with = {10^, ll^, 20^}. For the ID lat- 
tices we used kT = {0.5, 1.0, 8.0}; for the 2D lattices we used kT = {2.0, 2.1, ...3.0}, 
which includes temperatures on either side of T^. We also simulated a wider range of 
temperatures for = 100. We used the non-linear least squares method to fit the MC 
data to the model, allowing Agg and A to vary freely. Each data point was weighted by 
The errors in Agfi and A are 90% confidence intervals. For sufficiently high temper- 
atures, the model gives a very good fit to the data (see figures 4 and 5). In fact, some 
of the values are significantly lower than expected, which is likely due to equation 7 
giving an overestimate of the errors when F is close to 1. However, as the temperature is 

+ This models option (1) in the case of a 50/50 split. 
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decreased, the fit becomes poorer. This occurs at T ^ 0.3 for D = 1 and at T ~ 2.1 for 
D = 2, although the transition is not abrupt. At temperatures close to this transition, 
A^eff < 1 was often produced by the fitting, which is clearly unphysical. Additionally, 
the Gaussian approximation of the binomial distribution breaks down for the low values 
of Aeflf needed at low T. In general, however, the fitting process gave 1 < A'efr < A^ 
and < A < 0.5 as expected. The deviation of the actual lattice evolution from our 
model at low temperatures is not unexpected, as the model is extrapolated from the 
non-interacting case. However, we can see from figure 5(d) that a simple exponential 
model is also inadequate. 

We immediately notice that for low T the decay of F is significantly slower for 
D = 2 than for D = 1, when A^ and T are the same in each case (see figures 4 and 5). 
By considering the variation of the model parameters with A^, a deeper insight into the 
differences between the ID and 2D lattices may be achieved. (Variation of the model 
parameters with T is discussed in the appendix.) 



Fidelity of 1 D lattice with N=1 00 and kT=2.5 
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Figure 4. MC Data (blue with vertical error bars) and model (red) for ID lattices 
with (a) N = 100, kT ^ 2.5, A = 0.1707 ± 0.0003, iVcff = 46.7 ± 0.3, = 1093, 
Degrees of Freedom = 1334 and (b) N = 100, fcT = 8, A = 0.4005 ± 0.0011, 
iVefF = 99.9 ± 1.2, = 1421, Degrees of Freedom = 1570 
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Fidelity of 2D lattice with N=100 and kT=2.5 
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Fidelity of 2D lattice with N=100 and kT=2 
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Figure 5. MC Data (blue with vertical error bars) and model (or exponential fit for 
(d) only) (red) for 2D lattices with (a) = 100, kT = 2.5, A = (1.521 ± 0.010) x IQ-^, 
TVeff = 5.37 ± 0.09, = 71, Degrees of Freedom = 116; (b )N ^ 100, kT = 8, 
A = 0.282 ± 0.002, N^s = 62.2 ± 1.9, = 198, Degrees of Freedom = 198; (c) 
N = 100, fcT = 2, A = (8.52 ± 0.19) x 10-^ N^g = 0.80 ± 0.02, = 684, 
Degrees of Freedom = 358 and (d) N = 100, fcT = 2, A (1.113 ± 0.012) x lO'^, 
x'^ = 434850, Degrees of Freedom = 359. 



Figures 6a and 7a illustrate one of the key results of this work. For a ID lattice 
(figure 6a) A is independent of A^. For a 2D lattice (figure 7a), by contrast, A decreases 
as is increased. This is a manifestation of the fundamental difference in the dynamics 
of the ID and 2D lattices. In the ID lattice, no energy input is needed to increase the 
size of a region of incorrect spin, and hence increasing does not affect the rate of 
decay. In the 2D lattice, the opposite is true. Figure 7a also shows that the rate of 
decrease of A in the 2D lattice falls as we increase T. Our data are consistent with an 
exponential decay of A with increasing for T < T^. Interestingly, and in contrast with 
the ID case, for T > Tc there does seem to be some reduction in A as A^ is increased. 
However, figure 7a indicates that this reduction is subexponential in A^ for T > T^. 
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>. vs N for the 1 D lattice 
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N „ vs N for the 1 D lattice with l<T = 3.5 
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Figure 6. Variation in (a) A (at T = 2.0,4.0,8.0) and (b) N^e (at T = 3.5) with 
N for tlie ID lattice. Note some error bars are too small to be visible. We see that 
A is independent of N in this case, (b) shows the fitting Ncs = mN (red) with 
m = 0.61 ± 0.02. (c) shows the variation of m with T. 
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A vs N for the 2D lattice 



kT = 2.1 
kT = 2.2 
kT = 2.3 
kT = 3.0 



(a) 



N „ vs N for ttie 2D lattice with kT = 2.2 



I I 



So 100 150 200 250 300 350 400 450 
N 



(b) 



25r 
20- 
2^15 

10- 



N „ vs N tor the 2D lattice with kT = 2.7 



50 100 150 200 250 300 350 400 450 
N 



(c) 



Figure 7. Variation of (a) A (at T = 2.1,2.2,2.3,3.0, see dashed line to guide the 
eye) and (b),(c) A^cfi {T = 2.2 and T = 2.7) with N for the 2D lattice. These data 
are consistent with an exponential decrease in A as iV increases for T < Tc, and a 
sub-exponential decrease for T > Tc. 



For a ID lattice, Nes is approximately linearly proportional to A^: N^s — niN (see 
figure 6b). As T — )■ oo, m — )■ 1 as expected (see figure 6c). Therefore the scaling of the 
ID lattice with is not qualitatively different to that of A^ non-interacting spins: the 
only difference is in the value of m. The ID lattice maps directly to the non-interacting 
system as T — )■ oo. 

There are two competing factors influencing the variation of N^s as A^ is increased 
in the 2D lattice. Firstly, a larger lattice can accommodate more spin regions of a given 
size, which would lead to a positive correlation between A^efr and A^. Secondly, a larger 
lattice has more inherent stability due to interactions between spins, leading to to a 
negative correlation between A^efr and A^. Our simulations show that, just below the 
critical temperature, and A^ are approximately negatively correlated (flgure 7b), 
whereas, for T > Tc, Nea and A^ are approximately positively correlated (flgure 7c). 
Thus the phase transition appears to emerge naturally from the proposed model in the 
form of the variation of both A and A'efr with A^. We also notice that these correlations 
are not very strong for low T, indicating that the model is breaking down here. This is 
likely due to both the Gaussian approximation becoming poor, and the model becoming 
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4. Conclusions and Further Work 

Our model offers an intuitive view of the evolution of the Ising lattice at high and 
intermediate temperatures and provides continuity between the interacting and non- 
interacting cases: as we take the temperature to infinity, the model parameters approach 
their non-interacting values. The work described here entails a proof of principle of the 
proposed model, and an exploration of what we can infer about the ID and 2D Ising 
lattices from its parameters. Our model gives a quantitative backing to the commonly 
held belief that at low temperatures the 2D Ising lattice is a good physical memory, 
whereas the ID Ising lattice is not. More specifically, in the 2D case A decreases as we 
increase the system size, giving us significant returns on our investment in extra spins. 
In the ID case, by contrast, A is independent of system size so we get very little return 
on an investment in extra spins. Furthermore, we have seen the phase transition of the 
2D lattice emerge from our model, as the decrease in A as AT increases appears to change 
from exponential to sub-exponential as we move from T <Tc to T > Tc. 

Both classical and quantum data storage are likely to be of theoretical and practical 
interest in future years, prompting the study of the fidelity of a range of physical 
memories. This article begins this process with the study of the Ising lattice, the 
archetypal classical physical memory. 

A natural extension of this work would consider the applicability of the new model 
for /i 7^ 0, for D > 2, and for different boundary conditions. The low temperature 
case, for which neither the model proposed here nor an exponential model is entirely 
satisfactory, also merits further study. 
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Appendix A. Variation of Model Parameters with T 

The variation in the model parameters with T are qualitatively similar for the ID and 
2D lattices. As T is increased, A — )> 0.5 and A'cfr — )■ A^ as expected (see figure Al). The 
primary advantage of the 2D lattice is therefore in the scaling of the fidelity with A^. 
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Figure Al. Variation in (a),(c) A and (b),(d) N^ff with kT for the ID and 2D lattices 
with N = 100. Note some error bars are too smah to be visible. 
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